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We present a systematic calculation of the coherent multiple parton scattering with several nucleons 
in lepton-nucleus and hadron-nucleus collisions. We show that in £ + A reactions coherence leads to 
nuclear shadowing in the structure functions and a modification of the QCD sum rules. In p + A 
reactions we evaluate the nuclear suppression of single and double inclusive hadron production at 
moderate transverse momenta. We demonstrate that both spectra and dihadron correlations in p + A 
collisions at RHIC are sensitive measures of such dynamical nuclear attenuation effects. 



1 Introduction 

Copious experimental data 1 from central 
Au + Au reactions at the Relativistic Heavy 
Ion Collider (RHIC) has generated tremen- 
dous excitement by pointing to the possible 
creation of a deconfined state of QCD with 
energy density as high as 100 times normal 
nuclear matter density 2 . In order to better 
understand the jet quenching mechanism in 
A + A collisions we first need to address the 
multiple scattering between a partonic probe 
and the medium in simpler strongly interact- 
ing systems, like I + A and p + A. 

In this talk we present a systematic cal- 
culation of the coherent multiple parton scat- 
tering with several nucleons in lepton-nucleus 
and hadron-nucleus collisions. Recent the- 
oretical developments 3 ' 4,5 have been able to 
provide a consistent picture of the dynami- 
cal nuclear shadowing of sea quarks, valence 
quarks and gluons. Coherence leads to small 
x < 0.1 suppression 6 of the structure func- 
tions Ff 4 , and in I + A reactions 3,4 
and suppression 7 of single and double inclu- 
sive hadron production cross sections at small 
and moderate transverse momenta in p + A 
reactions 5 . Coherent multiple scattering and 




Figure 1. From 5 . Coherent multiple scattering of the 
struck parton in deeply inelastic scattering (a) and in 
hadron-nucleus collisions (b). 

the well-studied elastic multiple scattering 8 
co-exist in nuclear collisions and their relative 
role depends on the probes (or observables) 
and the underlying dynamics. The interplay 
between these two scattering channels is be- 
yond the scope of this talk. 

2 Coherence and Shadowing 

Hard scattering in nuclear collisions requires 
one large momentum transfer Q ~ xP 3> 
Aqcd with parton momentum fraction x 
and beam momentum P. A simple exam- 
ple, shown in Fig. 1(a), is the lepton-nucleus 
deeply inelastic scattering (DIS). The effec- 
tive longitudinal interaction length probed by 
the virtual photon of momentum is charac- 
terized by 1/xP. If the momentum fraction 
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Figure 2. From 3 . All-twist resummed F2(A)/F^(D) 
calculation versus DIS data 6 on nuclei. Data-Theory 
Ae/_t is also shown. 



of an active parton x <C x c = l/2rnjv^o ~ 
0.1 with nucleon mass tojv and radius ro, it 
could cover several Lorentz contracted nucle- 
ons of longitudinal size ~ 2ro(w7v / P) in a 
large nucleus. In the photon-nucleus frame, 
Fig. 1(a), the scattered parton of momen- 
tum £ will interact coherently with the par- 
tons from different nucleons at the same im- 
pact parameter 3 . In terms of the collinear 
factorization approach 9 , we showed 3 that the 
hard partonic interactions for any number of 
coherent multiple scattering in DIS are in- 
frared safe. The on-shell condition for £ fixes 
the incoming parton's momentum fraction to 
the Bjorken variable x = xb = Q 2 /(2P ■ q) 
but additional final state scattering forces the 
rescaling 3 
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with the scale of power corrections £ 2 sw 
lim^o \xG{x, Q 2 ). For £ 2 = 0.09 - 
0.12 GeV 2 the calculated reduction in the 
DIS structure functions, known as nuclear 
shadowing, is consistent with the xb~, Q 2 - 
and A-dependence of the data 6 as demon- 
strated in Fig. 2. 

Neutrino-nucleus DIS is a particularly in- 
structive example since it provides the pos- 
sibility to study the interplay between the 
dynamical power corrections and the heavy 
quark mass and is a unique handle on va- 
lence quark shadowing. For v + s — > £~ + c, 
the rescaling of the Bjorken-x reads: 
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where M c is the charm quark mass. From 
Eq. (2) it is evident that the role of the co- 
herent multiple final state scattering is to 
generate a dynamical parton mass mjjLj, — 
£ 2 (AV 3 _ i) j n th e nuclear chromomagnetic 
field 10 . The rescaling of the value of x B re- 
flects the conversion of the parton's energy 
into mass and correspondingly reduces the 
partonic flux. This attenuation depends on 
the x— slope of the parton distributions and is 
different for sea quarks and valence quarks . 

3 Suppression in Hadron-Nucleus 
Collisions 

Unlike in DIS, all diagrams with either final- 
state and/or initial-state multiple interac- 
tions in hadron- nucleus collisions, shown in 
Fig. 1(b), could in principle lead to medium 
size enhanced power corrections. However, 
once we fix the momentum fractions x a and 
zi, the effective interaction region is deter- 
mined by the momentum exchange q^ L = 
(x a P a — P c /zi)^. In the head-on frame of 
q — Pb, the scattered parton of momentum 
£ interacts coherently with partons from dif- 
ferent nucleons at the same impact parame- 
ter. Interactions that have taken place be- 
tween the partons from the nucleus and the 
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Figure 3. From 5 . Top: suppression of the single and 
double inclusive hadron production rates in d + Au 
reactions versus pj- for different rapidities. Bottom: 
the calculated nuclear modification at different im- 
pact parameters. The trigger hadron pyi = f -5 GeV, 
3/1 = 3 and the associated hadron y2 = 0. 



incoming parton of momentum x a P a and/or 
the outgoing parton of momentum P c / z\ at 
a different impact parameter are much less 
coherent and actually dominated by the in- 
dependent elastic scattering 8 . 

Similarly to the DIS case 3 , we find 5 that 
resumming the coherent scattering with mul- 
tiple nucleons is equivalent to a shift of the 
momentum fraction of the active parton from 
the nucleus in Fig. 1(b), 
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with the hard scale t = q 2 = (x a P a — P c /zi) 2 
and the color factor Cd depending on the 
flavor of parton d. Cq(q) = 1 and C g = 
Ca/Cf = 9/4 for quark (antiquark) and 
gluon, respectively, since the gluons couple 
twice as strongly to the medium via their 
average squared color charge. The shift in 
Eq. (3) leads to a net suppression of the cross 
sections, and the i-dependence of this shift in- 
dicates that the attenuation increases in the 
forward rapidity region. 

Dynamical nuclear effects in multi- 
particle production can be studied through 
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Centrality dependence is implicit in Eq. (4) 
and the modified cross section per average 
collision da^ hn /{N^ u ) can be calculated 5 
from the nucleon-nucleon cross section with 
the shift in Eq. (3). We consider d + Au re- 
actions at RHIC at VS = 200 GeV in the 
following discussions. 

The top panels of Fig. 3 show the rapid- 
ity and transverse momentum dependence of 
R dAu(b) and RdAuiP) for minimum bias col- 
lisions. The amplification of the suppression 
effect at forward y\ comes from the steep- 
ening of the parton distribution functions at 
small Xb and the decrease of the Mandelstam 
variable (— t). At high pti,PT2 the attenu- 
ation is found to disappears in accord with 
the QCD factorization theorems 9 . The bot- 
tom panels of Fig. 3 show the growth of the 
nuclear attenuation effect with centrality. 

Dihadron correlations ^(A^) = 

W~- — dA^ ass0Cia ted with 2^2 partonic 
hard scattering processes, after subtracting 
the bulk many-body collision background, 
can be approximated by near-side and away- 
side Gaussians. The acoplanarity, Aip ^ it, 
arises from high order QCD corrections and 
in the presence of nucler matter - transverse 
momentum diffusion 8 . If the strength of the 
away-side correlation function in elementary 
N+N collisions is normalized to unity, dy- 
namical quark and gluon shadowing in p + A 
reactions will be manifest in the attenuation 



of the 



Atfar — R 



hih- 



(by- 



iFar — ^ <- p A 

The left panels of Fig. 4 show that for 
Pt x — 4 GeV, pt 2 = 2 GeV the dominant 
effect in ^(At^) is a small increase of the 
broadening with centrality, compatible with 
the PHENDC 11 and STAR 12 measurements. 
Even at forward rapidity, such as y\ =2, the 
effect of power corrections in this transverse 
momentum range is not very significant. At 
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Figure 4. From 5 . Centrality dependence of C2(A</?) 
at various rapidities and moderate (left) and small 
(right) transverse momenta. Central d+Au and p+p 
data from STAR 12 . 

small =1-5 GeV, pr 2 — 1 GeV, shown 
in the right hand side of Fig. 4, the appar- 
ent width of the away-side (^(A^) is larger. 
In going from midrapidity, 2/1=2/2 = 0, to 
forward rapidity, y\ = 4, 2/2 = 0, we find a 
significant reduction by a factor of 3 - 4 in 
the strength of dihadron correlations. Pre- 
liminary STAR results 13 are consistent with 
our predictions. 

4 Conclusions 

In conclusion, we presented a systematic 
approach 3, 4:5 to the calculation of coherent 
QCD multiple scattering and resummed the 
nuclear enhanced power corrections to the 
structure functions F^, and F% and to 
the single and double inclusive hadron pro- 
duction cross sections in p + A reactions. At 
low Q 2 or — t oc p 2 -. we find a sizable sup- 
pression, which grows with centrality and via 
—t and Xb with rapidity. At high Q 2 ot —t 
the nuclear modification disappears in accord 
with the QCD factorization theorems 9 . We 
demonstrated that both particle spectra and 
dihadron correlations are sensitive measures 
of such dynamical attenuation of the parton 
interaction rates in p + A reactions. 

Our approach, with its intuitive and 
transparent results, can be easily applied to 



study the nuclear modification of other phys- 
ical observables in p + A reactions. The sys- 
tematic incorporation of coherent power cor- 
rections provides a novel tool to address the 
most interesting transition region between 
"hard" and "soft" physics in hadron-nucleus 
collisions. 
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